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ABSTRACT 

To investigate tlie poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel- 
ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500/im sources brighter than 50 mjy 
with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches 
at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron 
peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower 
frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. 
One of the two blazars is identified with a Fermi/LAT 7-ray source and a WMAP source. The physical parameters of the two 
blazars are briefly discussed. These observations demonstrate that the H-ATLAS survey will provide key information about the 
physics of blazars and their contribution to sub-mm counts. 

Key words. BL Lacertae objects: general - quasars: general - Submillimeter 



] 1. Introduction 

, Blazars, comprising BL Lac objects and flat-spectrum radio 

■ quasars, are a class of active galactic nuclei characterized by 
' remarkable properties: high luminosities at the two extreme 
, ends of the electromagnetic spectrum (i.e., in the radio and 
' 7-ray bands), rapid variability, apparent superluminal jet 
I speeds, flat or inverted radio spectrum, and high polariza- 
. tion in the optical waveband. They are interpreted in terms 
' of relativistic effects due to jets propagating at velocities 
, close to the speed of light almost along the line-of-sight. 

■ Their spectral energy distribution (SED) can be accounted 
' for as a combination of synchrotron emission, peaking at 
, a frequency i/^ that varies from ~ 10^^ Hz to > 10^^ Hz 
' (Nieppola et al. 2006), and inverse Compton scattering, as- 
, sumed to be responsible for the 7-ray peak. 

Padovani & Giommi (1995) classified BL Lac objects 
into 3 subclasses (low, intermediate, and high frequency 
synchrotron peaked BL Lacs) depending on the value of Vp . 
Abdo et al. (2009) extended the classification to all blazars, 
subdivided into low (LSP; i^p < 10^^ Hz), intermediate 
(ISP; 10" <iyp< 10^5 Hz), and high (HSP; z/^ > lOi^Hz) 
synchrotron peaked blazars. For many years, the blazar se- 
lection has relied mostly on either radio or X-ray surveys. 



* Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consor- 
tia and with important participation from NASA. 
** e-mail: gnuevo@sissa.it 



Radio surveys preferentially select LSPs, while X-ray sur- 
veys favor the HSPs. The resulting global distribution of 
Vp is thus bimodal, although objects peaking at intermedi- 
ate frequencies have also been found. The high sensitivity 
and nearly uniform sky coverage of the Fermi Gamma-ray 
Space Telescope made it a powerful tool for providing large 
blazar samples (Abdo et al. 2010). 

In the present paper, we report on the blazars selected 
at sub-mm wavelengths in the 16 deg Herschel- P^hAB 
(H-ATLAS; Eales et al. 2010) science demonstration phase 
(SDP) field. The Herschel data cover the poorly explored 
frequency range close to the synchrotron peak of the most 
luminous LSPs (Fossati et al. 1998; Padovani et al. 2006); 
Up carries crucial information about key physical parame- 
ters, such as the Lorentz factor 7^ of emitting electrons, 
the Doppler factor 5, and the magnetic field strength B 

Evolutionary models (De Zotti et al. 2005) predict 
~ 0.15 blazars deg~ brighter than 50mJy at 500 ^m, the 
approximate 5ct detection limit of the H-ATLAS survey 
(Rigby et al. 2010, in preparation), so that the full sur- 
vey may be expected to yield a sample of ~ 80 blazars 
selected, for the first time, at sub-mm wavelengths. This 
prediction assumes a flat radio spectral index a {a ~ 0.1; 
Si, oc z^~") followed by a parabolic decline close to and 
above a synchrotron peak frequency that increases with 
decreasing radio luminosity (Fossati et al. 2008). However, 
spectral measurements have provided evidence of a substan- 
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Table 1. The two H- ATLAS blazars found by cross-matching with the FIRST catalog. Coordinates (RA and Dec) refer 
to the H-ATLAS source. A (arcsec) is the angular separation between FIRST and Herschel positions. The quoted errors 
in Herschel fluxes are statistical only. Calibration errors are estimated to be ~15% for SPIRE (Swinyard et al. 2010), 
~10% at 100 /im, and ~20% at 160 /im. Upper limits are at the 3a level. 



H-ATLAS 
name 


A 5'l.4GHz 

arcsec nijy 


•S'sOOfim 

mJy 


•SaSOMin 

mJy 


S25011U1 
mJy 


•S'lTOMin 

mJy 


•S'lOOMin 

mJy 


J090910.14-012135" 
J090940.3+015957'' 


0.3 571.7 
10.3 317.5 


265.8 ± 8.9 
68.2 ±9.1 


193.8 ±7.2 
55.0 ±7.2 


159.6 ±6.4 
35.5 ±6.4 


105.4 ±26.5 
<76.2 


137.3±26.5 
<79.5 


" fHB89l 0906+015, 


"PKS 0907-f 022 
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Fig. 1. SED of [HB89] 0906±015 (non-simuhaneous data), 
in the observer's frame. Black dots are data taken from the 
NASA/IPAC Extragalactic Database (NED), red diamonds 
refer to H-ATLAS data, red circles are the Fermi/LAT data, 
red asterisks are the new ATCA data, red squares are the 
new Spitzer data, red stars are UKIDSS data, red hexa- 
grams are the GMRT data, and green triangles are 3cr upper 
limits derived from IRAS maps. The blue solid/dashed lines 
represent the best fit to the synchrotron/inverse Compton 
part of the SED, omitting the optical/UV fluxes, which are 
probably dominated by the direct emission from the accre- 
tion disk (see text). 



tial spectral curvature downwards from cm to mm wave- 
lengths (Waldram et al. 2007; Sadler et al. 2008). Vieira 
et al. (2009) find that their highest significance 150 GHz 
selected radio sources are consistent with fiat spectral be- 
havior (or a ~ 0.1) between 5 GHz and 150 GHz but a 
possible steepening between 150 and 220 GHz. Voss et al. 
(2006) reported an excess (based however on the detection 
of only three blazars) above extrapolations from lower fre- 
quency counts of the surface density at S* > lOmJy at 
1.2 mm, suggesting that at least a fraction of blazars have 
inverted (i.e. rising with frequency) spectra at mm wave- 
lengths, perhaps due to variability. The H-ATLAS data will 
also shed light on this issue. 

In this paper, we describe (Scct.[2|) the identification of 
blazars among sources detected in the H-ATLAS SDP field 
and discuss (Sect. ^ the properties of the 2 blazars found. 
In Sect. 2] we summarize our main conclusions. 



2. Blazar candidates 

Observations were performed with the ESA Herschel Space 
Observatory (Pilbratt et al. 2010). The 4deg x 4deg H- 
ATLAS SDP field, centered on (09''05'"30", +00°30'00"), 
was observed with the Spectral and Photometric Imaging 
Receiver (SPIRE; Griffin et al. 2010) at 250, 350, and 
500 /xm and with the Photodetector Array Camera and 
Spectrometer (PACS; Poglitsch et al. 2010) at 100 and 
IGOyitm. Details of the SPIRE and PACS mapmaking are 
given by Pascale et al. (2010, in preparation) and Ibar et 
al. (2010, in preparation), respectively. Sources were ex- 
tracted from the 250 /im map (see Rigby et al. 2010, in 
preparation). For each source, 350 and 500 /im fiux densi- 
ties were estimated by extracting fiux from the appropriate 
noise- weighted beam-convolved map at the source position 
determined in the 250 /im map. Fluxes at 100 and 160 /im 
were then assigned by matching to > 3ct PACS sources 
within a positional tolerance of 10". 

The selection of blazar candidates among H-ATLAS 
sources brighter than 50 mJy was made by identifying 
FIRST counterparts within a search radius of 15", the 
(Th ~ FWHM/2.35 of a Gaussian model of the 500//m beam. 
The probability that a true counterpart has an apparent 
positional offset > A is 



p(>A)=exph0.5(A/a)2], 



(1) 



where a is the global positional uncertainty, taking into 
account both H-ATLAS and FIRST positional errors. 
Sources detected with signal-to-noise ratio (S/N) greater 
than 10 at 250 /im, such as the blazar [HB89] 0906+015 
(J090910.1±012135; see Table[I]), have a cth-atlas - 1-4". 
For PKS 0907+022 (J090940.3+015957), which has data of 
lower S/N at 250 /xm, ctr-atlas is ^ 4.9". The FIRST po- 
sitions are accurate to better than 1 arcsec at the survey 
limit and to better than 0.5 arcsec for sources brighter than 
5 mJy at 1.4 GHz (White et al. 1997). With our choice of 
the search radius, the probabilities of having missed true 
associations at separations A > 15" are ~ 8 x 10"^'^ and 
~ 10"^, respectively, for the two values of cth- atlas and 
fpiRST — 0.5". We note that the two values of cth-atlas 
are representative of sources with data of high and low S/N 
at 250 /im. 

Since the distribution of FIRST sources is to a large 
degree uniform, the probability that one of them lies by 
chance within the angular radius A from a given Herschel 
source is 



= 7rA^ni.4GHz(> S), 



(2) 



where S is the flux of the sources we match to. The 
FIRST catalog is complete down to 0.75 mJy/beam. The 
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corresponding surface density is ni.4GHz(> 0.75 mJy) ~ 
89 deg~ , implying that the maximum probability of a 
chance association within 15 arcsec is 5 x 10~^. Since 
there are 185 sources with S'soopm > 50mJy in the H- 
ATLAS SD field, the expected number of chance associ- 
ations is ~ 0.9. Simulations made by offsetting the ra- 
dio positions at random and re-running the cross-matching 
yielded ric — 0.9(+2.4, —0.8) (68% confidence errors), con- 
sistent with expectations for Poisson statistics. 

The cross-matching yielded 19 matches. Of them, 10 
have far-IR/sub-mm colours typical of low-z dusty galax- 
ies. The far-IR/sub-mm SEDs of 4 others are consistent 
with those of dusty galaxies at z ~ 1-2 and their radio 
emission obeys the radio-FIR correlation (Yun et al. 2001). 
The other three sources are more than 4cr (a being the 
combined rms positional uncertainty) away from their near- 
est FIRST source, and are therefore very unlikely associa- 
tions. We note that having 3 spurious associations within 
our search radius is fully consistent with the results of our 
simulations. All of these 17 sources were discarded. We note 
that we aim to select only sources with non-thermal emis- 
sion brighter than 50 mJy at 500 /im, so we are not inter- 
ested here in the weaker blazar nuclei that may be hiding in 
dusty galaxies. The FIRST and Herschel photometric data 
for the 2 remaining sources are shown in Table [TJ 

Both H-ATLAS J090910.1-h012135 ([HB89] 0906-f015) 
and H-ATLAS J090940.3-I-015957 (PKS 0907+022) are 
known blazars. The former is a flat-spectrum quasar with 
strong broad emission lines at a measured redshift of 1.018 
(Falomo et al. 1994). The latter is a BL Lac with an esti- 
mated photometric redshift of 1.575 (Richard et al. 2004); 
a Keck spectrum obtained by E. Barton and J. Cooke 
does not show any clearly identifiable line. The surface 
densities of 1.4 GHz sources brighter than these sources 
are ~ O.ldeg"^ and ~ ldeg~^, respectively. According to 
Eq. (21), the probability that they fall by chance within 
their angular distance. A, from the H-ATLAS source is 
~ 2 X 10~^ and ~ 2 x 10~^, respectively, while, after Eq. ([T]), 
the probability that the true counterparts to the H-ATLAS 
sources have an apparent angular separation as large as the 
measured one, or larger, is 99% and 12%, respectively. 



3. Spectral energy distributions of [HB89] 
0906+015 and PKS 0907+022 

The photometric data for the two previously known blazars 
are shown in Figs. [T] and [2] To the data available from 
NED, we added the H-ATLAS fluxes, the UKIDSS data 
(Lawrence et al. 2007), new ATCA measurements around 
20 GHz, new GMRT measurements at 332 MHz (Jarvis 
et al., in preparation), new Spitzer data (Jarvis et al., in 
preparation) for [HB89] 0906+015 along with 3a upper lim- 
its that we obtained from IRAS maps (see Tabled. The 
ATCA measurements were performed by M. Massardi and 
L. Bonavera, in the framework of the Planck- ATCA Coeval 
Observation (PAGO) project. 

The blazar [HB89] 0906+015 was detected by the 
WMAP satellite (source 215 in the NEWPS5yr catalog; 
Massardi et al. 2009) and also by the Fermi/Large Area 
Telescope (LAT) (Abdo et al. 2010; note that in the pre- 
vious bright gamma-Ray source list the 7-ray source was 
identified with PKS 0907+022). The shape of its SED ex- 
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Fig. 2. SED of PKS 0907+022 (non simultaneous data) in 
the observer's frame. Black dots: data from the NED; red 
diamonds and triangles: H-ATLAS flux measurements and 
3(7 upper-limits, respectively; red asterisks: ATCA mea- 
surements; red five-pointed starts: UKIDSS data; red six- 
pointed starts: GMRT data; green triangles: IRAS 3(7 up- 
per limits. The blue solid line is the best-fit relation for the 
synchrotron part of the SED. 

Table 2. GMRT, ATCA, UKIDSS, and Spitzer fiux den- 
sity measurements and IRAS 3(7 upper limits for the two 
blazars. 



Flux (mJy) 


[HB89] 0906+015 


PKS 0907+022 


S332MHZ 


941 


144 


ftsGHz 


1101 + 55 


97 + 5 


5'24GHz 


742 ± 37 


69 + 4 


<S'l.03/jm 


0.38 


0.22 


Ol.25/im 


0.97 


0.24 


5*1.63^111 


0.47 


0.34 


5*2.20^111 


0.80 


0.54 


53.6(im 


1.46 + 0.09 


- 


'^4.5^im 


2.05 + 0.10 


- 


55.8/im 


2.80 + 0.12 


- 


Sg.OMm 


4.2 + 0.2 


- 


524/^m 


12.1+0.3 


- 


5l2/^m 


< 23 


<38 


•^25^111 


<42 


< 58 


•^60^111 


<40 


< 38 


5l00fim 


< 84 


< 145 



hibits the typical double hump distribution, attributable to 
synchrotron and inverse Compton, respectively. 

To determine the SEDs, we followed the approach of 
Abdo et al. (2009). First we fitted the part of the SEDs 



Table 3. SED parameters for the two blazars. Peak fluxes. 



i/_F(i/), in ergcm ^s ^. 



Name 


[HB89J0906+015 


PKS0907+022 


logio('^p)[HzJ 


12.6 


13.3 


\o^^^(ylF{vl)) 


-11.5 


-12.0 


logio(/^p^^)[Hz 


20.9 


- 


\o^^^,(ylpF{vi'')) 


-10.3 


- 
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dominated by the synchrotron emission using a third de- 
gree polynomial. From these fits, we obtained an estimate 
of the synchrotron peak energy, v^^ and peak intensity, 
i^^Fif^). Since [HB89] 0906+015 has strong broad emis- 
sion lines, its optical-UV flux is most likely dominated by 
the thermal emission produced by the accretion disk, with 
only a minor contribution from the beamed non-thermal 
continuum (Ghisellini et al. 2010b). There are indeed clear 
indications of an optical-UV bump. To avoid biasing esti- 
mates of the synchrotron peak frequency and intensity, we 
excluded these data from the fit. Abdo et al. (2009) demon- 
strated that the peak of the inverse Compton SED com- 
ponent, I'p^ , is strongly correlated with the 7-ray photon 
index (F). We used their best fit t'f'^-F relationship (their 
Eq. (5)) to estimate 1^^^ for the HB89 blazar (F=2.74). 
Finally we obtained the peak intensity by fitting the X-ray 
to 7-ray data points with a parabola peaking at u^'-' in the 
log(i/Fy)-logi/ plane. 

The SED parameters are summarized in Table [3l For 
PKS 0907-1-022, since no X-ray nor 7-ray measurements are 
available, we could only use the data fitting method for the 
synchrotron peak. In both cases, the available data from dif- 
ferent instruments are non-simultaneous. As for all blazars, 
they may therefore by affected by variability, which causes 
additional uncertainty in the derived spectral parameters. 
This adds to the uncertainties in estimates of spectral pa- 
rameters. Based on the best-fit model values of i^p, both 
blazars can be classified as LSPs; however, an ISP classifi- 
cation cannot be ruled out. 



wards from cm to mm wavelengths (Waldram et al. 2007; 
Sadler et al. 2008). 

The blazar catalog of Massaro et al. (2009) lists three 
more blazars in the SDP field: PKS0858-004 with S'i.4ghz = 
325 mJy, IRXS J085920.6-H0047 with S'i.4ghz = 39mJy, 
and IRXS J085749.8-h0135 with S'i.4ghz = 89mJy. None 
of them has a plausible counterpart among the H- ATLAS 
sources. In the case of PKS0858-004, this implies that the 
effective 1.4 GHz to 500 fim spectral index has to be steeper 
than 0.3, i.e., steeper than the effective indices of the two 
detected blazars (0.13 and 0.25), but still within the range 
observed for these sources. 

Based on the best-fit model estimate of the syn- 
chrotron peak frequency, v^, both [IIB89] 0906-1-015 and 
PKS 0907-1-022 fall in the low frequency synchrotron peaked 
(LSP) category. This suggests that the H-ATLAS sample 
will set strong constraints on the abundance of blazars with 
low (< 10^"^ Hz) values of i^p, which, in the blazar sequence 
scenario, are those with more powerful jets, more luminous 
accretion disks and higher black hole masses (Ghisellini et 
al. 2010a). 
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4. Conclusions 

The Herschel detection of two blazars in the SDP area is 
interesting in a number of ways. First of all, it confirms 
that the H-ATLAS survey can play an important role in 
determining properties of blazars in a crucially important 
spectral region. With an area of 550 deg , the full H-ATLAS 
should yield a sample of ~ 80 blazars, the first statisti- 
cally significant sample selected at sub-mm wavelengths. 
The large area is of crucial importance in this context, be- 
cause the blazar counts are expected to be relatively flat. 
The H-ATLAS survey will extend by about one order of 
magnitude downwards in flux the blazar number counts 
from the Planck survey, which is predicted to have a (95% 
reliability) detection limit of 430 mJy at 545 GHz (Leach 
et al. 2008). 

The number of detected blazars is consistent with the 
predictions of the De Zotti et al. (2005) model which yields 
~ 0.15 blazars deg~^ brighter than 50mJy at 500 /im, i.e., 
2.4 objects in the 16 deg~ SDP area. An observational es- 
timate of the blazar surface density at 5 GHz was obtained 
by Padovani et al. (2007; see their Figs. 4 and 6). Taking 
into account the steepening at frequencies approaching the 
synchrotron peak, the effective 5 GHz to 500 /im spectral in- 
dex yielded by the model is ~ 0.2. The 500 fim flux limit of 
50 mJy then corresponds to 130 mJy at 5 GHz; our results 
are then in very good agreement with the corresponding 
blazar surface density found by Padovani et al. (2007) and 
the relatively flat spectral index assumed by the de Zotti et 
al. model. If this result is confirmed by data of the full H- 
ATLAS blazar sample, future models will need to explain in 
addition evidence of a substantial spectral curvature down- 
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